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o : 

( We investigate the phase structure of dense QCD using the Schwinger-Dyson equation 

(SDE) with the improved ladder approximation in the Landau gauge. We use the gluon 
propagator with the electric mode corrected by Debye screening and the magnetic mode 
corrected by Landau damping. We solve the coupled SDE for the Majorana masses of the 
quark and antiquark (separately from the SDE for the Dirac mass) in the low density and 
intermediate density regions. In the low density region, both the SDEs for the Majorana 
masses and the Dirac mass have nontrivial solutions that correspond to the color symmetry 
breaking (CSB) vacuum and the chiral symmetry breaking (xSB) vacuum. Comparing the 
values of the effective potential in the two vacua, we show that the phase transition from the 
XSB vacuum to the CSB vacuum is of first order. The resultant value of the critical chemical 
potential is about 210 MeV, which is smaller than that obtained from the SDE analysis 
including only the Dirac mass. The momentum dependences of the Majorana masses of the 
quark and antiquark are shown to be quite different, although the masses are of the same 
order. 
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§1. Introduction 



Dynamical chiral symmetry breaking is one of the most important features of 
QCD. The QCD Lagrangian of the light quark sector has an approximate chiral 
symmetry, and it is spontaneously broken by the strong interaction of QCD. Several 
physical processes of light hadrons in the low-energy region are governed by the 
chiral symmetry property. In hot and/or dense matter, on the other hand, the chiral 
condensate (quark- antiquark condensate) melts, and the chiral symmetry is restored 
(for recent reviews, see, e.g., Ref. 1)). 

In early days of the research of such systems, it was pointed 2 )' 3 ) that there exists 
a diquark condensate (color superconducting condensate) in the high density region. 
Recently, it was shown 4 ) by using the instanton induced four- Fermi interaction model 
that the diquark condensate is on the order of 100 MeV. The diquark condensate is 
not a color singlet, and thus breaks the color gauge symmetry. This phenomenon is 
called color superconductivity (for recent reviews, see, e.g., Ref. 5)). 

There are many works on color superconductivity employing analysis based on 
the SDE (see, e.g., Refs. 6)-9)). Since the SDE analysis is valid for weak coupling, 
most of the works concentrate on the high density region. However, the density inside 
a neutron star is believed to be close to the phase transition point (see, e.g., Ref. 5)), 
where the competition between the chiral condensate and the diquark condensate 
becomes important. 
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Competition between the chiral condensate and the diquark condensate has been 
studied with several approaches (see, e.g., Refs. 10), 11), 9)). In Ref. 9) the SDE is 
converted into an algebraic equation using the confining model gluon propagator, 
and it was shown that the color symmetry breaking vacuum exists as a solution of 
the SDE, even in the low density region, but it is less stable than the chiral symmetry 
breaking vacuum. It is interesting to study whether such a false vacuum exists by 
fully solving the SDE. 

In most analyses previously carried out with the SDE, only one Majorana mass is 
included: The Majorana mass for the quark is analyzed while that for the antiquark 
is dropped, or the Majorana masses of the two are assumed to be equal. Thus, it 
would also be interesting to include two Majorana masses in solving the SDE. 

In this paper, we study the phase transition in QCD in dense matter using the 
SDE with the improved ladder approximation in the Landau gauge. In a dense 
medium, although the electric gluons have Debye masses, the magnetic gluons are 
not screened. The unscreened magnetic gluons give a dominant contribution to the 
formation of the Majorana mass gap. 12 -*' 6 -*' 7 -* In the high density region, through 
the hard dense loop approximation, the effect of the Debye screening is included 
in the electric mode of the gluon propagator, and the effect of Landau damping is 
in the magnetic mode, as shown in Ref. 6). Here we extrapolate the same form 
to the intermediate density region, and solve the coupled SDE for the Majorana 
masses of the quark and antiquark separately from the SDE for the Dirac mass. The 
true vacuum is determined by comparing the values of the effective potential at the 
solution. 

The symmetry breaking pattern in the color superconducting phase depends on 
the number of light quarks. When the number of massless quarks is two (Nf = 2), 
the dominant diquark condensate is the singlet of SU(2)l x SU(2)r, and thus the 
chiral symmetry is not broken, while the color SU(3) C is broken down to its subgroup 
SU(2) C (2SC phase). 4 ) For three massless quarks (Nf = 3), on the other hand, the 
dominant diquark condensate breaks the color S77(3) c and the chiral SU(3)l x 
SU(3)r symmetry together into their diagonal subgroup S , [/(3) c _|_£+/j [color flavor 
locking(CFL) phase). 13 -* The breaking pattern in real QCD depends on the strange 
quark mass (see, e.g., Ref. 14)). The above-mentioned two breaking patterns are 
understood as two limits of the strange quark mass: When the strange quark is 
heavy enough, the 2SC phase is realized. In contrast for the massless strange quark, 
the CFL phase is realized. 

In the present analysis we regard u and d quarks as massless, but we consider 
the current mass of the s-quark to be large enough to neglect the strange quark in 
the formation of the diquark condensate of u and d quarks; we assume that the color 
superconductivity is realized in 2SC phase. Furthermore, we ignore the effect of the 
strange quark in the effective potential, since we are interested in the differences 
among the values of the effective potential corresponding to the CSB vacuum, the 
xSB vacuum and the trivial vacuum. In the running coupling, however, the effect of 
the strange quark is included, because its mass is not larger than the energy region 
in which the dominant part of the support of the effective potential lies. In other 
words, we include the strange quark only quark in the present analysis. 
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The Lagrangian of the u and d quark sector in the present analysis is invariant 
under the chiral SU(2)l x SU(2)r symmetry. In the low density region this chi- 
ral symmetry is broken down to its diagonal subgroup, but the color symmetry is 
unbroken: 

SU(3) C x SU{2) L x SU(2) R -► SU(3) C x SU(2) V . 

In the high density region, on the other hand, the chiral symmetry is restored but 
the color symmetry is broken down to its subgroup: 4 ) 

SU(3) C x SU(2) L x SU(2) R -► SU(2) C x SU(2) L x SU(2) R . 

This paper is organized as follows. In §2, we summarize the quark propagator, 
the gluon propagator and the running coupling that we use in the present analysis. 
Several approximations of the quark propagator are made. We also give formulas for 
calculating the diquark condensate and the chiral condensate. In §3, we present 
the effective potential for the quark propagator and then derive the Schwinger- 
Dyson equation as a stationary condition of the effective potential. Section 4 is 
the main part of this paper, where we give the results of the numerical analysis of 
the Schwinger-Dyson equation. Finally, we give a summary and discussion in §5. In 
the appendices we summarize several intricate expressions and useful formulas. 

§2. Preliminaries 

In this section we briefly present the quark propagator, the gluon propagator and 
the running coupling that we use in the numerical analysis. In §2.1, we introduce the 
eight-component Majorana spinor (Nambu-Gorkov field) and give the general form 
of the full quark propagator as a matrix in the Nambu-Gorkov space. We obtain 
constraints on the full propagator from the parity and time reversal invariances. We 
further make several assumptions to restrict the form of the propagator. The gluon 
propagator with screening mass effects is presented in §2.2. We also give the explicit 
form of the running coupling that we use in our numerical analysis. We give formulas 
to calculate the quark-antiquark condensate and the diquark condensate in §2.3. 

2.1. Nambu-Gorkov fields and quark propagator 

Since we are interested in the phase structure of QCD, including the color su- 
perconducting phase, it is convenient to use the eight-component Majorana spinor 
(Nambu-Gorkov field) instead of the four-component Dirac spinors. The Nambu- 
Gorkov field is expressed as 

*=M*°)- * c=c * T - (2,i) 

where C = ij 2 j°. (We suppress the color and flavor indices until later in the paper.) 
Using the Nambu-Gorkov basis, the inverse of the free quark propagator is given by 

!sr(p) ^(P0 + ,h°-^- (m _ m) » __ /). (, 2) 
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Similarly, the full quark propagator is written in matrix form in the Nambu-Gorkov 
space as 

The four components of the above full quark propagator are not independent. From 
the relation tp c = Cip T it can be shown that they satisfy 

Sfii(p) =C[S F22 (-p)] T C- 1 , (2-4) 
S F21 (p) = C[S F21 (-p)] T C- 1 , S F12 (p) = C[S F12 (-p)] T C~ l , (2-5) 

where indicates the transposed matrix. Furthermore, we easily find the following 
relations among the components of the inverse propagator: 3 ) 

i{S F \p)}u + i{S F \p)}22 = 7°[i{S F \p)}ii + i{S F \p)} 22 h° , (2-6) 
i{S F Hp)hi = l [i{S F \p)}i2h° ■ (2-7) 

Here t indicates the hermitian conjugate matrix. 

The QCD Lagrangian at finite density is not invariant under charge conjugation 
(C) due to the existence of the chemical potential fj,, while it is invariant under time 
reversal (T) and parity (V). Time reversal of the full quark propagator is given by 

SF(p°,p)^f[S F (p°,-p)] T f , f= ( 1 ° J , T = z 7 y 7 , (2-8) 

and the parity transformation is expressed as 

s F {p ,p)-^r°s F {p ,-p)r\ r°=^ ° _°^ ) . (2.9) 

For [x = 0, charge conjugation is not broken explicitly. The charge conjugation 
transformation is expressed as 

S F{ p)—>C[S F {-p)] T C-\ (5=(q (2-10) 

Let us obtain constrants to the spinor, color and flavor structures of the propagator. *) 
According to the spinor structure, the inverse full quark propagator is generally 
parametrized by eight bases : 

{ 1, 75, 7o, 7570, 7-P, 757 ■ P, 7o7 ■ P, 75707" P }■ (2T1) 

Let {Sp 1 (p)}u (I, J = 1,2) denote the I- J Nambu-Gorkov component of the in- 
verse full propagator.**^ By using the above bases, general forms of {Sp 1 (p)}n and 

When the chiral symmetry is not broken, the following discussion is essentially the same as 
that in Ref. 15). 

**) Note that in our notation {S 1 ^ 1 }/,/ / {Sfij} -1 - 
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{Sp 1 (p)}i2 can be expressed as 

i{S F \p)}n = -B(p) + A(p)( Po + M ) 7o + C(p)7 • p + £>(p) 7 o7 ■ P 

+Z\+ 5) (p)7 5 /l+ + ^ 5) (p)75^ + ~+(;p)757o^ + H ( - ) (p) 75 7o^ , 

(2-12) 

^^(p)}^ = ^ + (p) 7 5^p + ^-(p)75^ + ^ + (p)7570^p + S"(p)757o^ 

+■8(5) (p) + %) (p)7o + C(5) (p)7 • P + %) (p)7o7 • P , (2-13) 

where A~ and yl+ are the projection operators for the quark and antiquark in the 
massless limit: 

It should be noted that the above functions B, A, ■ ■ etc., do not carry spinor 
indices, while they are still matrices in the color and flavor spaces. We assume 
here that 0(3) symmetry (spatial rotation) is not broken, so that B, A, etc., 
are functions of po and p = \p\. To avoid notational complexity we write B{jp) in 
place of B(po,p). In the following discussion, we often use B(—p), which represents 
B(-po,p)- 

Equations (2-4) and (2-7) relate {Sp 1 (p)}22 and {Sp 1 (p)}2i to {Sp 1 (p)}n and 
{S'p 1 (p)}i2, respectively. By using the notation in Eqs. (2-12) and (2-13), {Sp 1 (p)}22 
and {Sp 1 (p)} 21 can be expressed as 

i{S F \p)}22 = C[i{Sp\-p)} ll ] T C- 1 

= -[B(-p)f + [A(-p)] T ( Po - /i) 7 o + [0(-p)] T 7 • p + [ J D(-p)] T 7 o7 ■ P 
+[^ 6) (-p)]%4J- + ^^ 5) (-p)] r 75^- 

+ [~+(-p)f 7s7o /l p - + [S ( - 5) (-p)]%7 /l+ , (2-15) 

H- s 'f 1 (p)}21 = 7o[i{S F 1 (p)}i2]ho 

= - [A~(p)]^ 5 A+ - [A + (p)]^ 5 A p + [^+(p)]t 757o yl+ + [~-{p)]^^A- 

+ [£(5)(p)] f + [%)(p)] f 7o + [C ( 5)(p)] t 7-P- [%)(p)] f 7o 7 - P • 

(2-16) 

In the above expressions, the full quark propagator includes sixteen functions, 
which are still matrices in the color and flavor spaces. Here we obtain the constraints 
on these functions from Eqs. (2-5) and (2-6). First, Eq. (2-5) leads to 

A±(p) = [^(-p)] T , S±(p) = [^(-p)] T , (2-17) 

and 

5(5) (P) = [B {5) (-P)] T , A [5) (p) = -[A (5) (-p)f , 

C (5) (p) = [C (5) (-p)] T , D {5) (p) = [D {5) (-p)] T . (2-18) 
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Second, the relation (2-6) leads to 

(po + n)A{p) + (p - fi)[A(-p)] T = [(po + n)A{p) + (po - /i)[A(-p)] T ]t , 

B(p) + [B(-p)f = [B(p) + [ J B(-p)] T ] t , C(p) + [C(-p)] T = [C(p) + [C(-p)] T ]t , 

iZ>(p) + i[D(-p)] T = [iD(p) + i[D(-p)] T ] jf , (2-19) 

and 

i4+(p) + i[A+ ) (-p)] T = [iA^(p) +i[A( 5) (-p)] T ]i , 

Sfo (P) + [H ( - 5) (-p)] T = [H+(p) + [H ( - 5) (-p)] T ]t . (2-20) 

Let us further constrain the forms of the functions using time reversal (7~) and 
parity ("P). In various analyses, 2 ) " 4 )' 16 )' 17 ) it has been shown that the most favorable 
condensate carries even parity, though it was pointed out in Ref. 18) that a parity 
violating condensate can be formed. In the present analysis, however, we use the 
SDE with a one-gluon exchange kernel, so that there exists U(1)a symmetry in the 
system. This U(1)a symmetry enables us to choose the vacuum with even parity. In 
such a case, the number of the functions is reduced to eight: 

B (5) (p) = A {5) (p) = C (5) (p) = £> (5) (p) = , (2-21) 

A %) (P) = ^(5) W = E t) W = H ( _ 5) W = • ( 2 ' 22 ) 

If T is not broken spontaneously*^ , we obtain the following constraints on the func- 
tions allowed by V invariance: 

[A( P )] T = A(p) , [B(p)] T = B{p) , [C{p)] T = C(p) , [D(p)] T = -D{p) , 

(2-23) 

Now, let us consider the flavor and color structures of A, B, C, D, A^ and S' ± . 
As we discussed in the Introduction, we consider the strange quark mass to be large 
enough to ignore the strange quark in the formation of the diquark condensate of 
u and d quarks. Then color superconductivity is realized in the 2SC phase, where 
the color SU (3) c symmetry is broken down to its subgroup SU(2) C . In the chiral 
symmetry broken (xSB) phase, on the other hand, the chiral SU(2)l x SU(2)r 
symmetry is broken down to its diagonal subgroup SU(2)y. Functions in the full 
quark propagator are written to reflect these symmetry breaking structures. Here we 
consider the color and flavor structures in the mixed phase, where the color symmetry 
breaking and the chiral symmetry breaking discussed above occur simultaneously: 
SU{3) C x SU{2) L x SU{2) R -► SU(2) C x SU{2) V . The pure 2SC phase [SU(2) L x 
SU(2)r is left unbroken] or pure xSB phase [SU(S) C is left unbroken] can be realized 
in a certain limit of certain functions. The pure 2SC phase is realized for B = D = 

*-* Note that the effective potential shown in Eq. (3-4) is invariant under T transformation. 
Then, ImZl - = lmA + = is actually a trivial solution to the imaginary part of the SDEs (3-8) and 
(3-9). 
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S' ± = 0, and the pure xSB phase for = = 0. Under this breaking pattern, 
the possible color and flavor structures of the functions are given by 

A(p), B(p), C(p), D{p) ~ S ij (S ab - 5 a3 5 b3 ) , 5 ij 5 a3 5 b3 , 

A ± (p), S ± (p)^e tJ e ab3 , 8 tJ e ab3 , (2-24) 

where (a, b) and (i,j) are the color and flavor indices, respectively, and we fix the 
direction of the color symmetry breaking in the 3-direction, without loss of generality. 
Under the transpose of color and flavor indices, the function D must be symmetric, 
as dictated by the symmetry breaking pattern [see Eq. (2-24)], but antisymmetric, 
as dictated by time reversal invariance [see Eq. (2-23)], and hence the function D 
has to vanish: 

D(p) = . (2-25) 

Next we make several approximations of the structure of the quark propagator. 
At T = n = 0, the SDE with the ladder approximation in the Landau gauge leads 
to the fact that the wave function renormalization of the quark yields l(A(p) = 
C(p) = 1). Actually, this is necessary for consistency with QCD in the ladder 
approximation (see, e.g., Ref. 19)). For T > and/or jjl > it is also necessary to 
keep A(p) = C(p) = 1 by adopting, e.g., nonlocal gauge fixing. 20 ) In the high density 
region, however, the deviation of the wave function renormalization from 1 have been 
shown to be small for any choice of the gauge parameter. 6 ) In addition, the relations 
A(p)^l and C(p)^l do not imply chiral symmetry breaking. Therefore, in this paper 
we ignore the deviation from 1 and regard A(p) = C(p) = 1 as an approximate 
solution for general fi, even though we keep the Landau gauge. Furthermore, we 
assume that the chirality even channel in the majorana masses is dominant and 
ignore the chirality odd parts: 

S+(p) = ~-(p) = . (2-26) 

In fact, in the high density region, where the chiral condensate vanishes, two massless 
quarks with the same helicity condense in the scalar (J=0) channel(see, Ref. 15)). 
Finally, since the color and flavor structure of the diquark condensate is given by 
eije ab3 , we ignore the terms proportional to Sij€ ab3 in the function A^: We fix the 
color and flavor structures of A ± as 

A ± (p) ~ e^" 3 . (2-27) 

From the above discussion, the inverse full quark propagator can be expressed 



as 



iS F (p)~ 



(Po + - 7 • P ~ B{P) A(p) \ 

A{p) (Po- -7-P- B(-p) J ' 

(2-28) 



8 



M. Harada and S. Takagi 



where*- 1 

A(p)° b = e ij e ab3 l5 [A+(p)A+ + A-{p)A~\ 

r ab3 



-eae 



1 s[A+{p)A- + A-(p)A+] , 



B(p)$ = B 1 (p)5 ij (5 ab - 5 a3 5 b3 ) + B 3 (p)S ij S a3 S b3 . (2-29) 

Now the full quark propagator includes four scalar functions, B\ and B% correspond- 
ing to the Dirac masses responsible for the chiral symmetry breaking and A + and A~ 
corresponding to the Majorana masses responsible for the color symmetry breaking. 
From Eq. (2-19), the Dirac masses B\ and S3 obey the constraint 

Im[B i (p)] = -Im[B i (-p)] . (i = l,3) (2-30) 

Equations (2-17) and (2-23) imply that the Majorana masses A + (p) and A~(p) are 
real and even functions of po : 

A ± (p) = A±(-p)=[A ± (p)]* . (2-31) 

By taking the inverse of the expression in Eq. (2-28), the full quark propagator 
is given by 6 ) 

/ R^{p) - {(po + Ai)7° ~ 7 • P— B{p)}^ l A(p)RZ 1 {p) 

{P) { -{{po-^)l°-l-p-B(-p)}- 1 A(p)R- + \p) 




where 

R+{P) = {(Po + M)7° - 7 ■ P ~ B(p)} - A(p){( Po - /i) 7 ° - 7 ■ p - B{-p)Y l A{p) , 
R-{p) = {(po - fih° ~ 7 ■ P ~ B(-p)} - A(p){(p + / u) 7 ° - 7 ■ p - B(p)}"M(p) . 

(2-33) 

To use Sp in the Schwinger-Dyson equation, we need explicit forms of R+ 1 , RZ 1 , 
and so on. We list them in Appendix A. 

2.2. Gluon propagator and the running coupling 

In dense medium, gluons generally acquire Debye masses. On the basis of the 
renormalization group equation 12 ) and the Schwinger-Dyson equation, 6 )' 7 ) in the 
high density region it has been shown that the long range interaction mediated by 
magnetic gluons gives the dominant contribution to the formation of the Majorana 
mass gap. The value of the Majorana mass gap obtained from these analyses is 
consistent with that derived using models based on the instanton induced interac- 
tion. 4 ) In addition to the Debye masses, five gluons have Meissner masses, due to 

f) Note that in Eqs. (2-13), (2-16), (2-17), (2-23) and (2-24), A+ and A~ are still matrices in 
the color and flavor spaces. In Eq. (2-29) the same notation is used to express the scalar functions 
after factoring the color and flavor indices. 
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the Meissner-Higgs effect, when the color SU(3) is broken to its subgroup SU(2). 
According to the result of the SDE analysis in the high density region, 6 ) these masses 
are of no importance in the SDE. 

In this paper we include the Debye screening and the Landau damping in the 
gluon propagator through the hard dense loop approximation. However, we ignore 
Meissner masses. Moreover, as explained in the previous subsection, we take the 
Landau gauge for the gluon propagator and assume A(p) = C{p) = 1 which is 
needed for consistency with QCD in the ladder approximation (see, e.g., Ref. 20)). 

The explicit form of the gluon propagator that we use in this paper is given by 6 ) 

D^ B (k) ee 5 AB D, u (k) 

= i5 AB \M_ (1) + iS AB 1 (2) 

\k\3 + 7TM D 2 \k 4 \/2 (fc 4 )2 + | fc |2 + 2M 2 f«" 

(2-34) 

where k± = —iko and Mp is the Debye mass of the gluon. In the hard dense loop 
approximation this is given by 6 -* 



Md = ■ (2-35) 



0$,(i = 1,2) are the polarization tensors defined by 

(« • k ? ™ ^(2) _ (« • k) 2 



(1) = p± , V"' pu (2) = pu (2 .o 6) 

(u • k) 2 - k 2 " u ' " (u • k) 2 - k 2 ' 1 ; 



where 



p j_ _ _ k^ky _ kpk v ^ k^u u + u^k u UpU v 2 (C) o 7 \ 

r^-9v> fc2 > i >- fc2 u . k + ( u . fc )2 fc • 

The Lorentz four-vector u M = (1,0,0,0) in the gluon propagator reflects the explicit 
breaking of Lorentz symmetry due to the existence of the chemical potential in the 
rest frame of the medium. 

In the SDE at zero density, it is important to use the running coupling, since the 
high-energy behavior of the mass function derived from the SDE with the running 
coupling is consistent with that derived from the operator product expansion. 19 ) One 
consistent way to include the effect of the running coupling is to use the improved 
ladder approximation, 21 ) in which the high-energy behavior of the running coupling 
is determined by the one-loop renormalization group equation derived in QCD and 
the low-energy behavior is suitably regularized. In the present analysis we consider 
low and medium density regions, so that we use the Higashijima-Miransky-type 21 )' 22 ) 
running coupling 

n (F) - 9HE) _ WqiiVe - 2N f ) 

as{E) - — - max( M/ ) ' (2 ' 38) 
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where 

t = ]n-^-, t f = (2-39) 

^qcd ^qcd 

with E being the energy scale, yl qc( j the characteristic scale of QCD,*) and Ef the 
infrared cutoff scale introduced to regularize the infrared singularity. **) As discussed 
in the Introduction, in the present analysis we assume that the mass of the strange 
quark is large enough to ignore the s-quark in the formation of the diquark conden- 
sation of u and d quarks. On the other hand, it is natural to assume that the current 
mass of the s-quark is smaller than yl qc d- In such a case, the effect from the s-quark 
should be included in the running coupling. ***) Thus we set Nf = 3 and N c = 3 in 
the running coupling (2-38). Note that the value of the running coupling in the in- 
frared region needs to be sufficiently large in order to involve the effect of dynamical 
symmetry breaking, i.e., tf has to be sufficiently small. Here we set if = 0.5, around 
which the various physical quantities for jjl = are very stable with respect to the 
changes in tf. 22 ) In the numerical analysis we investigated the tf dependence of the 
results. 

2.3. Condensates 

In this subsection we give formulas to calculate the chiral condensate and the 
diquark condensate. 

The chiral condensate is generally expressed as 

d 4 p 



(f2\^t(0)\f2) A = -J Ptv[S F11 ], (240) 



(27T) 

where \Q) is the ground state at nonzero density and the trace is taken in the spinor, 
flavor and color spaces. Summations over the color index a and the flavor index i are 
implicitly taken on the left-hand side of Eq. (2-40). The quantity A is the ultraviolet 
cutoff, introduced to regularize the logarithmic divergence. In the actual numerical 
analysis we introduced two cutoffs for the temporal and spatial components of the 
momentum, but they are expressed symbolically as A in this subsection. For the 
present form of the full quark propagator in Eq. (2-32), this is given by 

(n\j>Wm\n)A 
<- A d 4 p 



i(2vr) 4 
16 



F{p,B u A) 




Hf-tf ~ {Bi(- P )} 2 Bi(p) - A + (p)A-(p)B 1 (-p) 



Here, yl qcd is determined from the infrared structure in the present analysis, while the usual 
Aqcd is determined from the ultraviolet structure. 

**' When we study the high density region, we need to add another cutoff scale, ln(/i//l qc d), at 
which the running stops. However, this is neglected here, because we study only the intermediate 
region around the phase transition point. 

***' As we show in the next section, yl qcd = 604 Me V in the present analysis, which is obviously 
larger than the s-quark mass. 
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+ F( P , bIa = 0) ~ ~ p2 ~ W-P)} 2 ) *3(P)] , (2-41) 
where F is defined as 

F(p, B 1 ,A) = [(p + p) 2 - f - {B 1 (p)} 2 ][( Po - /i) 2 - f - {^(-p)} 2 ] 
-[(Po) 2 - (p-v) 2 }\A + (p)\ 2 - [(pa) 2 - (p + v) 2 ]\A-(p)\ 2 
+\A + (p)\ 2 \A'(p)\ 2 + 2B 1 (p)B 1 (-p)A+(p)A-(p) . (2-42) 

For the chiral condensate of the charge conjugated quarks, we obtain the following 
relation from Eq. (2-4): 



f A d4 P , [Q 1 f A d4 P t rc 



(2vr) 4 y (2tt) 

H^KY>f(0)l^ • (2-43) 

These expressions take the following familiar form when we set [i = 0, A = and 
Si = S3 = 5: 



4AT TVf f ^ (2-44) 

4iVciV/ i i ( 2 ^-(p r+pi + {B(p)y ' i244j 



^ d 4 p B(p) 

where iV c = 3 and N f = 2. 

The diquark condensate is generally expressed as 

{Q^ea^^C^m^A = ~f ^tr[ e ( c ) e (% F1275 ] , (2-45) 

where and e*^ are antisymmetric matrices in the color and flavor spaces, respec- 
tively: 

{e (c) } ab = e ab3 {e (/) } .. = e .. (2 . 46) 

In the present approximation, the diquark condensate is given by 
(^|(e^e a63 )^ T ]?C 7 5^(0)|^)^ 
d A p 8 



(Po) 2 - (P - - {^"(P)} 2 - |£i(p)| 2 U + (P) 



+ |(po) 2 - {P + ti? ~ {A + (p)} 2 - ^(prf^A'ip) 

(2-47) 

Using time reversal invariance [see Eq. (2-8)], we obtain 

{m^^M]tc^c] b j(o)\n) A 

= (^ tr[e(c)e(/) ^ 2l75] = S'^-M^^Sf^] 

= -(m T ]?C^ b A0)\[2) A . (2-48) 
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We note that the parity invariance existing in the present analysis leads to a vanishing 
parity violating condensate: 



{n\(eVe ab3 )[r]?Ctf(0)\f2) A 



/ /1 (2^ tr[e(C)e(/)5F12] = - (249) 



(27T) 



In the improved ladder approximation at zero density, the high-energy behavior 
of the mass function is consistent with that derived using the operator product 
expansion (OPE). The chiral condensate calculated using the mass function was 
shown to obey the renormalization group evolution derived with the OPE (see, e.g., 
Refs. 19)). Then, we identify the condensates, which are calculated with cutoff 
A, with those renormalized at the scale A in QCD. Therefore we scale them to 
the condensates at 1 GeV, using the leading renormalization group formulas. The 
relation between the chiral condensate at the scale A and that at the scale E is given 
by 



ln£/7l qcd 
In A/A qcd 



</2|#^(0)|rt>, 



where 



K 



9C 2 (F) 
llN r - 2N 



N 2 

c 



/ 



117V„ - 2N f 2N, 



(2-50) 



(2-51) 



Noting that the attractive force between two quarks in the 3 channel due to one-gluon 
exchange is one half of that between a quark and antiquark in the singlet channel, 
we immediately obtain the following relation between the diquark condensate at the 
scale A and that at the scale E: 



(n\(e l h ab3 M T ]?ctf(o)\n) 



In E/A qcd 
ln/l/yl qcd 



k/2 



(n\(e l h ab3 M T ]«cMo)\n) A 



(2-52) 



§3. Effective potential and the Schwinger-Dyson equation 

In this section we present the effective potential for the quark propagator and 
then derive the Schwinger-Dyson equation (SDE) as a stationary condition of the 
effective potential. 

As explained in the Introduction, we assume that the current mass of the strange 
quark is large enough to ignore it in the valence quark sector, and we include only 
u and d quarks in the effective action. Then, the effective action for the full quark 
propagator Sf is given by 23 ) 

r[S F ] = \ ^-iTrhn(S^) - ^S^Sf) - ir 2PI [S F ]j , (3-1) 

where Tr and Ln are taken in all the spaces and / 2pi[5i?] represents the contributions 
from the two-particle irreducible (with respect to the quark line) diagrams. The 
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factor of 1/2 appears because we are using the eight-component Nambu-Gorkov 
spinor basis. In the high density region, the one-gluon exchange approximation 
is valid, since the coupling is weak. In the present analysis, we extrapolate this 
approximation to intermediate density and include only the contribution from the 
one-gluon exchange diagram in I^pi^f]) 



1 



r 2P i[S F ] = ~-Tv(S F • igFl ■ S F • igl% ■ D$?) , (3-2) 
where is the quark-gluon vertex in the Nambu-Gorkov basis, defined as 

F " A = ( 7 A -Y{T A ) T ) ' (3 ' 3) 

From the effective action (3-1), the effective potential in the momentum space can 
be written 



v[s F ] = -r[s F }/ J d 



4 x 



^4 (in det{S F (p)} - tv{Sp-\p)S F (p)} 

4^ r A± 



-II 

(3-4) 

where In, det and tr are taken in the spinor, color and flavor spaces. 

The SDE is obtained as the stationary condition of the effective potential 
(5V[S F ]/5S F = 0): 

S- 1 = SP" 1 - (zgr^ ■ S F ■ igTl) • . (3-5) 

In principle, Eq. (3-5) leads to coupled equations for seven functions [B, A, C, A ± 
and S ± in Eqs. (2-12)-(2-16)] in the quark propagator, which are still matrices in the 
color and flavor spaces. As discussed in §2.1, we assume that A = C = 1 and E ± = 
are the approximate solutions of the equations. Furthermore, in the present pattern 
of the symmetry breaking, we take the color and flavor structures of the Dirac and 
Majonara masses as to be those given in Eq. (2-29). As a result, the above SDE (3-5) 
leads to four coupled equations for four scalar functions, B±, B%, A + and Z\~. These 
are given by 

Blip) = J -^^a s D, u ( q -pM^T A S FU (qh u T A 5{ c) ] , (3-6) 
B 3 (p) = J -^7Ta s D, v ( q -pM^T A S F11 (q)YT A 8i c) ] , (3-7) 
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^-(P) = / j^lna s D^(q-p)tT[rT A S F12 (q)Y(T A ) T A; l5 e^e^] , 

(3-8) 

A+{P) = I ^4^«^^(9-p)tr[ 7 ,i ^5 F12 (g)7^(r A ) T 4 75e ( c )6W] , 

(3-9) 

where 

{5[ c) } a6 = <5 ab - S a3 5 b3 , {4 C) } a6 = £ a3 ^ 3 , (3-10) 

and the traces are taken in the spinor, flavor and color spaces. The quantity a s = 
a s {E) on the right-hand sides of Eqs. (3-6)-(3-9) is the running coupling defined in 
Eq. (2-38). For consistency with the chiral symmetry, the argument of the running 
coupling should be chosen as E 2 = — (p — q) 2 - 20 - ) However, it has been shown 20 ) that 
the angular averaged form E 2 = —p 2 — q 2 is a good approximation at T = fx = 0. 
For this reason, in the present analysis we use the angular averaged form for \i > 0. 
Graphical representations of these gap equations are displayed in Fig. 1, and the 
explicit forms of these SDEs are given in Appendix A. 



B(p) = 




Fig. 1- Graphical representations of the gap equations. 

As we show in Eq. (2-31), the Majonara masses A + (p) and A~(p) are real 
and even functions of po, while in general the Dirac masses B\(p) and B%(p) are 
complex functions. Equation (2-30) is the constraint on the imaginary part, while 
no constraint is obtained on the real part from general considerations. However, as 
shown in Refs. 24) and 25) for the case of A + (p) = A~(p) = 0, the structure of 
the SDE leads to a natural constraint on the real parts of the Dirac masses B\{p) 
and B^{p). This is shown as follows. Let us take the complex conjugate of the 
SDEs for B\ and B3 in Eqs. (A-4) and (A-5). By using the fact that the kernel Kq 
and the Majorana masses A + and A~ are real functions, it is easily shown that the 
combination {Bf(p),B^(p)} satisfies the same coupled equations as {B\(p*), Bz(p*)}. 
Since p^* = (p°*,p) = {—ipi,p) = {—p°,P), this means that {B\(p), B^(p)} satisfies 
the same coupled equations as {Bi(—p), B^(—p)}. Thus these two combinations are 
equal up to their signs, {Bf(p),B|(p)} = {B^-p), B 3 (-p)} or {Bf(p),B|(p)} = 
— {Bi(— p), Bs(— p)}. Since the general relation (2-30) implies that imaginary parts 
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of B\{p) and B 3 (p) are odd functions of p°, the above argument implies that the real 
parts are even functions of p°: 



(3-11) 



Using the above properties for B\ and B 3 and those for A + and A~ in Eq. (2-31), we 
can always restrict the integration over the temporal component of the momentum, 
P4 = — ip°, to its positive region. 

Substituting the solution of Eq. (3-5) into Eq. (3-4), we obtain the effective 
potential at the vacuum, i.e., at the stationary point. Since the effective potential 
itself is divergent, we subtract the effective potential at the trivial vacuum: 

V^[2\ + ,2\-,Bi,B 3 ] 
= V[A + ,A-,B 1 ,B 3 ]-V[0, 0,0,0] 



r _d^p 

J i(2^ 



(2tt)' 



d 4 p 
i(2ir) 4 



2 In 



F(p,B 1 ,A) 



+ ln 



[ipo + M) 2 - P 2 ] [ipo ~ M) 2 - P 2 ] J 

' F(p,B 3 ,0) ~ 

[(Po + M) 2 - P 2 ][(Po - m) 2 - P 2 } 



F(p,B 1 ,A) 



[(po - m) 2 - P 2 ~ {Bi(-p)} 2 ][(po + - P 2 } 



+ 



1 



F(p,B 3 ,0) 



+[(P0 + V? -P 2 ~ {^i(p)} 2 ][(P0 - ti? -P 2 } 
-[(Po) 2 ~ (p + n) 2 ]\A~\ 2 - [(po) 2 - {p-n) 2 ]\A+\* 

( [{po ~ V) 2 -P 2 - {B 3 (-P)} 2 } [(po + V? - P 2 ] 

+ [(P0 + M) 2 " P 2 " {^3(P)} 2 ][(P0 " ^ " P 2 ] 



-6 



(3-12) 



Here F is defined in Eq. (2-42). The value of the effective potential in Eq. (3-12) 
is understood as the energy density of the solution. Therefore the true vacuum 
should be determined by evaluating the value of the effective potential. The vacuum 
corresponding to the solution becomes more stable as the value of V^ Q i decreases. 

§4. Numerical analysis 

In this section we give the results of our numerical analysis for 0.2 < [i/A qc< i < 
1.0. The parameters necessary to carry out the numerical analysis are ^l qc d and the 
infrared cutoff parameter tf in the running coupling. /l qc d is the unit of the energy 
scale in our numerical analysis, and it is determined by calculating the pion decay 
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constant f n for fixed tf at zero density through the Pagels-Stokar formula, 26 ^ 



/, = 4^;?e*e — biTew — ■ (4 ' 1) 

We use f w = 88 MeV in the chiral limit 27 ) as an input. For [i = 0, the dependences 
of the physical quantities on tf have been shown to be slight around tf = 0.5. 22 ) 
With this in mind, for the time being, we fix tf = 0.5 for general \x. From these 
inputs we obtain A qc< i = 604 MeV. Later in this section we study the tf dependence 
of our results. 

We introduce the framework of our numerical analysis in §4.1. Then, we give 
the solutions for the Majorana masses (A~ and A + ) fixing B\ = B3 = in §4.2. 
We give the solutions for the Dirac mass {B\ and -B3) with A~ = A + = in §4.3. 
Finally, in §4.4, we study the phase transition from the hadronic phase to the color 
superconducting phase. 

4.1. Framework of the numerical analysis 

In this subsection we summarize the framework of our numerical analysis. First, 
as discussed below Eq. (3-11), we note that it is always possible to restrict the p^- 
integral (p^ = —ipo) to its positive region due to the properties in Eqs. (2-31) and 
(3-11). To solve the SDEs numerically, we transform the variables p^ and p = \p\ 
into new variables U and X. For these transformations, we use density-independent 
transformations in the low density region (pL < /iq) and density-dependent transfor- 
mations in the medium density region (/j, > fio), where /j,q is determined in §4.2. 
In the low density region, where the chiral condensate is formed, the characteristic 
scale of the system is vl qc d- The dynamical information comes mainly from the region 
P4, p < /lqcd- On the other hand, in the high density region, where the diquark con- 
densate is formed, the chemical potential /x, in addition to ^l qc d> gives an important 
scale. Here, the dynamically important region is p^ < A qc d and p ~ fi. Therefore we 
adopt the transformations p^ — > U and p — > X determined by the following: 

Pa = A}cd • exp(£7) , p = yl qcd • exp(X) for fi < /i , (4-2) 

P4 = 3fi • exp([7) , p = 3/U • exp(X) for /j, > /j,q . (4-3) 

Under these transformations, the integrations over p$ and p on the interval [0, 00] 
are converted into integrals over U and X on the interval [—00, 00]. In the numerical 
integration, we introduced ultraviolet (UV) and infrared (IR) cutoffs for U and X, 
restricting their values as 

U€[A IR ,Auv], X G [Xir, Xuv] ■ (4-4) 

We discretize these intervals of U and X evenly into Njj and Nx points, respectively: 

U[I] = A IR + AU ■ I, J = 0,1,---,J%-1, 
X[J\ = Xi R + AX ■ J, J = 0,1,---,N X -1, 

(4-5) 
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where 

AU = AU A V ~ A ' R , AX = X ^~ X ! R . (4-6) 
Nu-1 N x -1 v ' 

The integrations over p4 and p are thus replaced with the following summations: 

[ d P4 ^AUj2e U[I \ [dp^AXj2e X[J] - (4-7) 
•* I J 

In the present analysis, for the UV and IR cutoffs, we use 

p A : {A IR , A uv ) = (-12.5, 2.5) , p : (X IR , \ uv ) = (-3.5, 2.5) . (4-8) 

The validity of these choices is checked below. 

We solved the SDEs with an iteration method. Starting from a set of trial 
functions, we updated the mass functions with the SDE: 



{5 li3old , Z\± d } => [Right-hand sides of SDE's (3.6)-(3.9) | => {Bi i3mw , ^L} • 

(4-9) 



Then we stopped the iteration when the convergence condition 
d A p 1 



eA% d > I T^TTtr 



sv \ ] ( sv 



5[S F (p)} \S[S F (p)] 



= J -0^ ^2\B lold (p) - B lncw (p)\ 2 + \B md (p) - B 3ncw (p)\ 2 

+2\A+ d ( P ) - Ai cw (p)\ 2 + 2\A- lA {p) - ZC w (p)| 2 j (4-10) 

was satisfied, with suitably small e. In the present analysis we set e = 1CP 10 . 
4.2. B = , A ^ solution 

In this subsection we report the results of the numerical solution of the SDEs 
with Dirac masses fixed to zero (B\ = B% = 0). The initial trial functions used here 
are 

Bi(p) = B 3 (p) = , 
A + ( P ) = 0, 

A-(p) = A qcd . (4-11) 

In this case, the outputs become B\ (p) , B% (p) = and A~(p),A + (p) / for all fj,. 
We call this solution the color symmetry breaking (CSB) solution. 

To check the validity of the UV and IR cutoffs in Eq. (4-8) we plot the solutions 
A~(p) and A + {p) and the integrand of V^oi at n/A^ = 0.70 in Figs. 2 and 3, 
respectively. Figure 2 shows that both A~ and A + become small in the UV region 
of p4, as well as that of p. It has been shown(e.g., in Ref. 6)) that A + is very small 
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compared with A~ in the high density region. In the medium density region, Fig. 2 
shows that A + is of the same order as A~ , but their momentum dependences are 
rather different. This implies that it is important to include two Majonara masses 
A + and A~ in the low density and intermediate density regions. From Fig. 3 we see 
that the dominant contribution to the effective potential lies within the integration 
range. These figures show that the choices of the ranges in Eq. (4-8) are sufficient 
at fi/Aqd = 0.70. We carried out similar analyses for all the cases we study in this 
paper, and confirmed that the choices of the ranges in Eq. (4-8) are sufficient for the 
present purposes. 




Fig. 2- Momentum dependences of Majorana mass functions A (p)/A„ cd 
and A + (p)/A qc d at fi/A qc d = 0.70. Here U = ln(p4/3/i) and X = 
ln(p/3//). 




Fig. 3- Integrand of V sol /A* cd at n/A qcd = 0.70. The upper 9/10 of the 
figure is clipped. Here U = ln(pi/3fi) and X = ln(p/3/x). 
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IVsoil 1/4 /Aqcd-sgn[V sol ] 



"30-30" + 

"40-40" x 

"50-50" * 

"60-60" □ 



|Ll/Aqcd 

Fig. 4- Typical values of (V^oil^Mqcd ■ sgn[V so i] for four choices of the 
fineness of the discretization, N X )=(30,S0), (40,40), (50,50) and 
(60,60). 



-<\JAJ/> 1/3 /Aqcd 



"30-30" + 

"40-40" x 

"50-50" * 

"60-60" □ 



Jl/A}cd 

Fig. 5- Typical values of — (V'VOiGev/^qcd f° r f° ur choices of the fineness 
of the discritization, (J%, N x )=(30, 30), (40,40), (50,50) and (60,60). 
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"30-30" 
"40-40" 
"50-50" 



0.3 0.5 0.7 0.9 M'/vVjcd 

Fig. 6- Typical values of A~(p4 = 0,p = /i)/M qc d for four choices of the 
fineness of the discretization, (N U ,N X )={S0,30), (40,40), (50,50) and 
(60,60). 



IVsoil 1/4 /Aqcd-sgn[v S oi] 



"schemel" 
"scheme2" 



H/A, 



cd 



Fig. 7- Typical values of iV^oil^Mqcd • sgn[V so i] in two different schemes 
of the discritization for 0.30 < /V^4qcd < 0.38. Here, "scheme 1" is that 
in which Eq. (4-2) is used and "scheme 2" is that in which Eq. (4-3) is 
used. 
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A~(ppo>P=M)/Aicd 

"schemer ~ 
"scheme2" x 



(l/Aqcd 

Fig. 8- Typical values of A~(p4 = 0, p = n)/A qc d in two different schemes 
of the descritization for 0.30 < fi/A qc( i < 0.38. 



Next, we elucidate the dependence of the results on the size of discretization. 

We show typical values of | V^soi 1 1/4 /^qcd ■ sgn[V so i]*) in Fig. 4, -(V'VOiGevAVd in 
Fig. 5 and A~(p^ = 0,p = [i)/A qc d in Fig. 6 for four choices of the fineness of the 
discretization, (Nu, N x )= (30,30), (40,40), (50,50) and (60,60). To obtain (^)iGeV, 
we used Eq. (2-52) with 



A 2 = A 2 qcd [exp(2A uv ) + exp(2X uv )] for n < /x , (4-12) 

A 2 = (3fi) 2 [exp(2A uv ) + exp(2X uv )] for ^ > Mo ■ (4-13) 



These figures show that the values Njj = 60 and Nx = 60 are sufficiently large for 
the present purpose. 

Let us next check that the two schemes of the discritization described by Eqs. (4-2) 
and (4-3) are smoothly connected with each other. For this purpose, we calculated 
the values of V^oi and Z\~(j>4 = 0,p = n) for several values of in 0.30 < [i/A^ < 
0.38 using the two schemes. We show the results in Fig. 7 and Fig. 8. These 
clearly show that the two schemes are smoothly connected with each other around 
MoMqcd = 1/3- For this reason, in the following analysis we fix ^o/^qcd = 1/3. 



*) sgn[V; o i] = 1 for V so i > and -1 for V sol < 0. 
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J/3, 



-<\|Ap 1/J /Aqcd 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



1/3 

Fig. 9- Chemical potential dependence of — (ipip) 1Ge y / A qc( \ 
for 0.2 < /X/Mqcd < 10. 



Jl/Aqcd 



Vsoil 1/4 /Aqcd- sgn[\^oi] 







CSB 



"Majorana" 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 MY Aqcd 

Fig. 10- Chemical potential dependence of iKoil^Mqcd ■ sgn[V so i] 
for 0.2 < fi/A qcd < 1.0. 
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Now, we display the resultant values of the diquark condensate in Fig. 9, where 
we have used Eq. (2-52) to obtain (ipi/i) iGeV- Figure 9 shows that (ipi/i) iGeV is 
between — (250 MeV) 3 and — (150 MeV) 3 , which is comparable to the value of the 
chiral condensate for fi = 0: (^V')iGcV = —(225 ± 25 MeV) 3 . 28 -* Furthermore, this 
figure shows that there exists a diquark condensate for all // we studied. 

As discussed in §3, the true vacuum is determined by evaluating the value of 
the effective potential at the solution. We plot the chemical potential dependence of 
iKsoil^VAjcd ■ s g n [V so i] in Fig. 10. This figure shows that the value of the effective 
potential is always negative, which implies that the CSB vacuum is always more 
stable than the trivial vacuum with B = A = 0. This is consistent with the result 
obtained in Ref. 9), where the SDE was converted into an algebraic equation by 
using the confining model gluon propagator. 



A(ppo>P=M)/^Vd 



"Gap" 




0.05 - 



Fig. 11- Chemical potential dependence of A (p4 = 0,p 
0.2 < p/A qcd < 1.0. 



JjVAqcd 

AOMqcd for 



Next, we plot the chemical potential dependence of A~(p^ = 0,p = ^)//l qc d 
in Fig. 11. This shows that the Majonara mass gap is on the order of 100 MeV 
in this low density region. We also plot the p dependence of A~(p4 = 0,p)/A qcd 
at /V^lqcd = O-^ an d 1-0 in Fig. 12 and the p^ dependence of A~(p4,p = /i)/Aj C d 
in Fig. 13. These figures show that the entire scale of the mass function becomes 
small when we increase /i. From Fig. 12, we see that A~ has a peak around p = /i 
(X = — In 3 = —1.1). However this peak is not sharp, and this implies that the 
effects of the quarks away from the Fermi surface are important in the medium 
density region. 
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A (p 4 =o)Mqcd 

0.15 



0.10 



0.05 




-ln3 
(_P=L0 

Fig. 12- p dependence of A (p)/M q cd at Pi = 
for r/ = /J/Mqcd = 0-5 and 1.0. Here X = ln(p/3//). 



X 



A (p=o)/Aqcd 

0.15- 



0.10 



0.05 



Air 



r(=0.5 



rj=1.0 



T| = [l/Aqcd 




-ln3 
(P 4 =^ 

Fig. 13- pa dependence of A~(p)/ A qc d at p = fx 
for r\ = /j/ylq C d = 0.5 and 1.0. Here U = ln(p4/3/x). 
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4.3. B ^ , A = solution 

In this subsection we report the results of the numerical solution of the SDEs 
with the Majorana masses fixed to zero (A^(p) = 0). The initial trial functions used 
here are 

Bi(p) = A qcd , 
B 3 (p) = , 

A + (p) = A~{p) =0 . (4-14) 

In this case, the outputs become B\(p) = B^(p), A^{p) = for all /i. We obtain 
Bi${p) 7^ for small fi, and B\-^{p) = for large [i. We call the former solution the 
chiral symmetry breaking (xSB) solution. 

In Fig. 14 we plot the solutions Re[l?i(p)] and Im[£?i(p)] (or equivalently Re[i?3(p)] 
and Im[i?3(p)]) at /X/Mqcd = 0.30. This figure shows that both the real part and the 
imaginary part become small above A qc & (U, X=0), and that the imaginary part 
has a peak around p4 = A qc d (U=0). These structures are consistent with those 
obtained in Ref. 25). 




Auv Auv 



Fig. 14- Momentum dependence of the Dirac mass function B\(p) / A qc d 
at fi/A qcA = 0.30. Here U = ln(p 4 /ylQCD) and X = \n(p/AqcT>)- 

We display the resultant values of — (V^)iGev/^q cd anc ^ |Kioi| 1//4 /^qcd - sgn[V so i] in 
Figs. 15 and 16, respectively. Considering Fig. 15 alone, it might seem that the phase 
transition occurs around /i/Mqcd = 0.40. (In Ref.24), the chiral phase transition point 
was determined by the point where their iteration converged to the trivial solution.) 
However, Fig. 16 shows that the value of V so \ is positive for 0.36 < /i/Aj C d < 0.40, 
although it is negative for ///Ajcd < 0.36. In other words, the value of the effective 
potential of the %SB vacuum for 0.36 < /VAjcd < 0.40 is larger than that of the 
symmetric vacuum. This implies that the xSB vacuum for 0.36 < /i/^4 qc d < 0.40 
is a false vacuum, and that the nontrivial solutions there correspond to metastable 
states. The existence of metastable states was found in Ref. 29) by assuming the 
momentum dependence of the mass function, and it was also shown in Ref. 25) by 
fully solving the SDE. The value of the critical chemical potential /x c /^qcd — 0.36 in 
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the present analysis is smaller than that obtained in Ref. 25), where the screening 
mass of the gluon was not included and a slightly different form of the running 
coupling was used. 

— 1/3 / A 
-<\|A|/> /Aqcd 





"Chiral" + 


++ XSB 








1 1 i 





0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



: Jl/Aqcd 

Fig. 15- Chemical potential dependence of (V^)iG C vMq C d ^ or ^-2 < 

/U/Mqcd < 1.0. 
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Fig. 16- Chemical potential dependence of l^oil^/Aicd ■ sgn[V so i] for 
0.2 < /i/A qcd < 1.0. 
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4.4. Phase transition from hadronic phase to superconducting phase 

In this subsection we compare the values of the effective potential for the CSB 
solution with that for the %SB solution to determine the true vacuum. We combine 
the effective potential for the CSB solution in Fig. 10 with that for the xSB solution 
in Fig. 16 and display them in Fig. 17. Figure 17 shows that, although the CSB 
vacuum is more stable than the trivial vacuum, the xSB vacuum is the most stable 
among these vacua in the low density region. Hence the true vacuum in the low den- 
sity region is the xSB vacuum. This is natural, because the strength of the attractive 
force between two quarks in the 3 channel is weaker than that between a quark and 
antiquark in the singlet channel. We find that the chiral phase transition and the 
color superconducting phase transition occur simultaneously at /J, c /A qc d = 0.344 and 
that the phase transition is of first order. We note that this phase transition occurs 
at lower density than the chiral phase transition in the absence of color supercon- 
ductivity. 

To compare the diquark condensate with the chiral condensate, we display them 
together in Fig. 18. Note that these are scaled to 1 GeV by the renormalization 
group formulas given in Eqs. (2-50) and (2-52). This figure shows that the diquark 
condensate for all n is of the same order as the chiral condensate in the low density 
region. The resultant values of these condensates at the critical chemical potential 
are (^V)iGeV = - (241 MeV) 3 and (^)iGeV = - (166 MeV) 3 . 



IVsoil 1/4 /Aqcd-sgn[v S oil 



"Dirac" 
"Majorana" 



x *xx Xx 



XXx *xx 



CSB 

xxx xxxx 



xxxxxxx: 



x xxxxxxxxxxxxxxxxxxxxxxxxxxx_ 

XSB 

Jl/A}cd 

Fig. 17- Values of |y so i| 1/4 Mqcd • sgn[V so i] (0.2 < n/A qcd < 1.0). Here 
"xSB" indicates that the inputs B / and A = were used, while 
"CSB" indicates that we used B = and A / as the initial trial 
functions. 
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-<\|n|/> 1/3 /Aqcd , -<\|n|/> 1/3 /Aqcd 



"Chiral" 
"Color" 



XSB 



CSB 



xxxxxxxx 



xxxxxxx 



xxxxxxxx 



xxxxxxx 



<xxxxxxxx 



0.9 

vV3 



|Ll/Aqcd 



Fig. 18- Values of -(^iGevMqcd and - {^)^ eV / A qcd (0.2 < 
/VAicd < 1-0). Here "xSB" indicates the chiral condensate and "CSB" 
indicates the diquark condensate. 

Finally, let us check the tf dependence of the critical chemical potential. We list 
the values of the critical chemical potential fi c for tf =0.3-0.7 in Table I, together 
with the values of the chiral condensate and the diquark condensate at fi c . This 
table shows that the value of \i c is almost independent of the value of tf. Although 
the value of (i/iip) depends slightly on that of tf, (ipip) is quite stable with respect to 
the changes in tf. 

Table I- Dependence of the critical chemical potential fi c on the infrared 
cutoff parameter tf. A qc d is the value determined from f n = 88 MeV 
through Eq. (4-1). The values of the chiral condensate and the diquark 
condensate shown together are at the critical chemical potential [i c . 



tf 


yl qcd [MeV] 


MMeV] 




-W>lSv[MeV] 


0.3 


568 


210 


223 


169 


0.4 


579 


207 


231 


168 


0.5 


604 


208 


241 


166 


0.6 


647 


214 


254 


167 


0.7 


712 


223 


269 


168 



§5. Summary and discussion 



We studied the phase structure in two- flavor dense QCD by solving the Schwinger- 
Dyson equations for the Dirac and Majorana masses of the quark propagator with 
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the improved ladder approximation in the Landau gauge. We showed that the mo- 
mentum dependences of the Majorana masses of the quark and antiquark are quite 
different, although the masses are of the same order. This implies that the inclusion 
of two masses is important for studying the intermediate density region. 

The values of the chiral condensate (quark-antiquark condensate) and the di- 
quark condensate were calculated using standard formulas. We identified the con- 
densates, which were calculated with the cutoffs A, with those renormalized at the 
scale A in QCD, and then scaled them to the condensates at 1 GeV using the leading 
renormalization group formulas. The resultant value of the diquark condensate is on 
the order of — (200 MeV) 3 , which is comparable to that of the chiral condensate for 
H = 0. 

The true vacuum was determined by comparing the values of the effective po- 
tential at the solution. We found that there exists a color symmetry breaking (CSB) 
solution for all the values of the chemical potential we studied, 0.2 < /x/Aicd < 1j 
and that the CSB vacuum is more stable than the trivial vacuum in all regions. 
Comparing the value of the effective potential for the CSB vacuum with that for the 
chiral symmetry breaking (xSB) vacuum, we showed that the %SB vacuum is more 
stable than the CSB vacuum in the low density region. We also found that the phase 
transition from the xSB vacuum to the CSB vacuum is of first order. The critical 
chemical potential was determined as \i c = 207-223 MeV. 

Finally, we make several comments. In an analysis employing the four-Fermi 
model, n ) it was apparently shown that there is a mixed phase, in which both the 
chiral condensate and the diquark condensate exist, in a small region of the chemical 
potential. In our analysis, we tried to find the solution corresponding to this mixed 
phase. However, we could not find such a solution in the present iteration method, 
although we used several initial trial functions. The value of the critical chemical 
potential obtained in the present analysis may depend on the detailed form of the 
running coupling. It is interesting to perform the similar analysis by using the 
different forms of the running coupling such as the one used in Ref.25). 

In the present analysis we included the effects of the Debye mass and the Lan- 
dau damping into the gluon propagator using the hard dense loop approximation, 6 ) 
assuming smooth extrapolation from the high density region. As we can see easily 
from the expression in Eq. (2-34), the magnetic mode explicitly breaks the Lorentz 
invariance, and therefore we cannot apply the form in Eq. (2-34) to the zero chem- 
ical potential case. It may be interesting to compare the present results with those 
obtained using the gluon propagator without effect of Landau Damping. This will 
be done elsewhere. 
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Appendix A 

Quark Propagator and Schwinger- Dyson Equation 



The explicit forms of the Nambu-Gorkov components of the full quark propaga- 
tor in Eq. (2-32) are given by 



Sfu(p)$ 



and 



Sf22(p)$\\m 



- M ;5(p) - B(-p);A ± (p) - -A*(p)] 



F(p,B 1 ,A) 

({po -P + A*)[(P0 " V? " " {^i(-p)} 2 ] -{Po-P- /u)|ZT (;p)| 2 ^ 70 
-Bi(-p)2\-(p)2\+(p) + B 1 {p)[(p - ^) 2 - P 2 - {Bi(-p)} 2 ]|yl+ 
+ | ( (po + P + /*o - /u) 2 -P 2 - {^i(-p)} 2 ] - (po +P- V)\A + ( P )\ 2 ho 



-Bi(-p)^+(p)2\-(p) +5i(p)[(po-/i) 2 -p 2 - {Bi(-p)} 2 ] \A: 



(po-p + n)-yo + B 3 (p) }A: 



5ijS 



ah 



{po+fl) 2_ p 2_ {B3(p)} 2 

+ \ (Po +P + IJ)lQ + B 3 {p) \A~ 



c%<5 a3 <5 63 



(A-l) 



Sfi2(p)% = S F2 i(p)%\[» - -m;S(p) - J B(-p);Z\ ± (p) - -^(p)] 
= -i{(po + ^7° - 7 • P ~ B(p)}- 1 A(p)R^(p)- 1 ^ 

= ^ ^ | ^(Po -p + /i)Z\ + (p)#i(-p) - (po -p- /i)^ - (p)Bi(p)^7o 
+ (\B l (p)\ 2 A + (p) - [(p ) 2 - (p + M ) 2 - \A + (p)\ 2 )A-(p)} y$ 

+ \{(po+P + l^)A-{p)B l {-p) - (p + p - i^)A + {p)B l (p) j 70 



+ (\B 1 {p)\ 2 A-{p) - [(po) 2 - (p- ^) 2 - |ZT(p)| 2 ]Z\+(p)) }^ 



afe3 



(A-2) 
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where 

F(p, B 1 ,A) = [(po + fi) 2 -P 2 - {B 1 (p)} 2 ][(p - /j,) 2 -p 2 - {B.i-p)} 2 ] 

-[(Po) 2 - (P-») 2 ]\A+( P )\ 2 - [(p ) 2 - (p + tf]\A-{p)\ 2 
+ \A+{p)\ 2 \A~(p)\ 2 + 2B 1 (p)B 1 (-p)A + (p)A-(p) . 

(A-3) 

Substituting the expression for Sfu in Eq. (A-l) into the SDEs for B\ and B% 
in Eqs. (3-6) and (3-7), we obtain 



(2vr) 4 
xKo(Q4,P4,q,P) 



dqq 2 a s 



B 3 (q) 



5 F + (q,B 1 ,A) 1 

6 F(q,Bi,A) + 2(iq A + ^-<?-{Bz{q)} 2 



(A-4) 



B 3 (p) = -2vr 



d<?4 / , . 

dqq a s 



oo (2vr) 4 



x-M<?4,P4,<?,p) 



F+(q,B 1 ,A) 1 flsfo) 

Ffai^) + 3(% + /x)2-^-{£ 3 (<7)} 2 



(A-5) 



where 

= #i (<?)[(% - P) 2 - 9 2 - {i?i(-9)} 2 ] - B 1 (-q)A + (q)A-(q) , 

(A-6) 

with ^4 = — igo and p4 = — ipo- The integration kernel i^o is given by 
K (q4,P4,q,P) =-iJ dQD^{q-p)ti{^A^Y) 



.- + p| 3 + -? +2 . (1 + d)4 r log k- + p| 2 +^ + 2M 2 



4 • — log F — _ |3 ^ 3 
Sqp \q — p^ + uf 



qp~'° \q-p\ 2 +lo 2 + 2M|, ' 

(A-7) 



where 



w ? = | M D^ i uj=\q A - p A \ . 



When we set A = A + = 0, we have 
F+(q,B 1 ,A = 0) 



Bi(q) 



F(q, B ± ,A = 0) (iq 4 + ~ q 2 ~ {B 1 (q)} 2 



(A-8) 



(A-9) 



If we further set B\ = B3, the two equations in Eqs. (A-4) and (A-5) become identi- 
fied, This implies that B\ = B3 is actually a solution of the SDEs for A~ = A + = 0. 
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Next, substituting Eq. (A-2) into the SDEs for A~ and A + in Eqs. (3-8) and 
(3-9), we obtain 



A 



■ 2 G + (g,B l ,A) _ 2 G_(q,B u A) 

Ki{q4,P4,q,P) ■ ~ ■ -pr, — 5 — tt- + K 2(q4,P4,q,P) ■ - 



3 F^B^A) 



3 F^B^A) 
(A-10) 



dqq 2 a s 



Ki{q4,PA,q,P) ■ o • EvI d — JV" +^2(94,1*4, 9, p) • o 



3 F{q,B 1 ,A) 



3 F(q,B 1 ,A) 
(A-ll) 

where 

G + (g, B l5 2\) = Bx^BxC-g)^^) + [(g 4 ) 2 + (g + + |^ + (g)| 2 ]ZT (g) , 

(A-12) 

G-(<?, -Bi,^) = ^(^(-^"(g) + [(g 4 ) 2 + (g- a*) 2 + |ZT (g)| 2 ]Z\+(g) • 

(A13) 

The integration kernels i^i and ii^ are given by 

Ki(qi,Pi,q,p) = ~ij df2D^(q - p)tr(A±^AjY) 
(q 2 - p 2 ) 2 + Z' min(g, p) 



= — 7T- 



y/3u?q 2 p 2 



arctan 



+ \q 2 — p 2 \ — ui max(g,p) 



-7T 



(q 2 -p 2 ) 2 - cof ui + |g + p| ( _ (g 2 - p 2 ) 2 - uf ^ aj + |g + p| 2 - o^jg + p| 



3ujfq 2 p 2 



uji + \q-p\ 



+ TT- 



ln 



6ujfq 2 p 2 cof + \q - p\ 2 - ui\q - p\ 



+7 r4lnj + lg ; + g! + * 



(g 2 -p 2 ) 2 (g + p) 2 
2M| ) + tc' 2 n (g-p) 2 



3gp + |g — p| 3 2g 2 p 2 
{(g + P) 2 + 2M 2 + w 2 }{(2M 2 + uj 2 ) 2 + (g - pfuo 2 } (g + p) 2 + 2M 2 + u? 



2M^{2Ml + 



In 



(g - p) 2 + 2M 2 + cj 2 



{(g+p) 2 +cj 2 }{(g - p) 2 +o; 2 } (g + p) 2 + uj 2 



2M 2 D 



(g -p) 2 +^ 2 



+ird 



_2_ _ (g - p) 2 + cj 2 (g + p) 2 + cj 2 
gp 2g 2 p 2 (q — p) 2 + uj 2 



(A-14) 



K 2 (q4,P4,q,p) = ~ij dQD^q - p)ti(A±>f Af-y") 
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(q 2 - p 2 ) 2 + uj 



ir- 



\f?>ujfq 2 p 2 



arctan 



^/3uJl mm(q,p) 



uf + 



oji max(g,p) 



(q 2 -p 2 ) 2 -uj? oji + \q+p\ _ ^ {q 2 - p 2 ) 2 - uf uf + \q + p\ 2 - u t \q + p\ 

uj? + \q-p\ 2 -ui\q-p\ 



3ujfq 2 p 2 
4 ,,3 



ln ^i i+lQ + P 



ui + \q-p 

3 



— 7T 



Qujfq 2 p 2 



In 



(q 2 -p 2 ) 2 jq+pf 
2M 2 +w 2 n (q-p) 2 



3qp u)f + \q— p\ 3 2q 2 p 2 

{(q - p) 2 + 2M 2 D + u> 2 }{(2M 2 D + uj 2 ) 2 + {q + p) 2 uj 2 } (q + p) 2 + 2M 2 D + uj 2 
" — > - " (q p) 2 ■ 2M 2 - J 2 



2M 2 



2Ml(2M% + uj 2 
In 



{(q + p) 2 + UJ 2 }{(q - p) 2 + uj 2 } (q + p) 2 +uJ 2 



-nd 



2 

qp 



{q+pf 



UJ 



2q 2 p 2 



In 



{q — p) 2 + uj 2 

{q + p) 2 + uj 2 
(q -p) 2 +uj 2 



We note that for B\ = we have 
G±(q,B 1 =0,A) 



-AT(q) 



F(q, B x = 0, A) (q,) 2 + (q T m) 2 + \^(q)\ 2 



(A-15) 



(A-16) 



Then the SDEs (A-10) and (A-ll) are identical to the well known forms given 
in, e.g., Ref. 6). 



Appendix B 

Convenient Formulas 



In this appendix we give several convenient formulas to obtain the explicit forms 
of the Schwinger-Dyson equations and the effective potential given in §3 and Ap- 
pendix A. 

As shown in §2.2, we use the Landau gauge in the present analysis, and two 
polarization tensors are used in Eq. (2-34). When we use the general covariant gauge, 
as done in Ref. 6), we need the three independent polarization tensors defined as 



(1) = p± + ( U ■ fc ) 2 pu (2) = (" • k ? pu (3) = kjJt* 

^ ^ (u ■ fc) 2 - fc 2 ^ ' » v " (u ■ fc) 2 - fc 2 ^ ' " fc 2 



where 



p _L _ k^k u 



The following formulas are convenient to obtain the SDEs: 

0«tr(yM±r) = 4 , Og)tr( 7 %V) = 2 , (i = 2, 3) 



(B-l) 



fc 2 ' ^~ fc 2 u-fc + («-fc) 2 • 1 j 
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0«tr(4 7 %V)= 2(l + t) 

o£M4^ A W = 2(1 - 1) 
oj?M4^ A h u ) = 2 (i - o 

0$tr(7l± 7 %^)=2(l + t) 



g 2 +p 2 - gp(l + t) 

g2 _|_ ^2 _ 2gpi 

g 2 + p 2 + gp(l-£) 

^2 _|_ ^2 _ 2gj)t 

g 2 +p 2 + gp(l - t) 
q 2 + p 2 — 2qpt 

q 2 +p 2 - qp{l + t) 
q 2 + p 2 — 2qpt 



~ (1 " t) 
~ (1 + 1) 



(q + p) 2 + (g4 -pa) 2 

q 2 +p 2 — 2qpt + (c/4 — P4) 2 

(g ~ P? + (g4 ~ Pi? 
q 2 +p 2 - 2qpt + (g 4 - p 4 ) 2 



O^tr^^T") = (l + t) 



g 2 + p 2 - 2gpi + (g 4 - p 4 ) 2 

(g-ff) 2 + (g 4 ~p 4 ) 2 

' g 2 + p 2 - 2qpt + (g4 - Pa) 2 



(B-3) 



where t = p • g/|p||g|. Angle integrations in the SDEs are performed using the 
following formulas: 6 -* 



r dL 

1 



dfl\q— p\ 



— l g^ + ^ |3+ ^ 3 



\q-p\ 3 +ujf 3qp |g-p| 3 + 



3 ' 



dfl 



7r . |g + p| 2 +cj 2 +2M 
log 1 



|g-pl 2 +w 2 + 2M 2 qp 6 |g-p| 2 +uj 2 +2M 2 D ' 



7T 



2 (g 2 -p 2 ) 2 +^ 
gp \/3u>fq 2 p 2 



arctan 



/ V^3aj min(g, p) 

\Lvf + |g 2 - p 2 \ - loi max(g, p) 



| (g 2 - p 2 ) 2 - wf + \q + p\ (g 2 - p 2 ) 2 - ujf ^ u 2 + |g + p| 2 -t^|g + p| 



3ujfq 2 p 2 



loi + \q-p\ 



I 4 ln ^f + k+P| 3 



3gp u;f + |g-p| 3 



Qujfq 2 p 2 



ujf + |g -p| 2 - w;|g -p| 



d.!?|g — p| 



= 7T 



_2_ _ (g 2 -p 2 ) 2 +w; 
gp 



0«tr(4 7 'M3V) 



\/3wi min(g,p) 



\[?>J 2 q 2 p 2 



arctan I .-1 

V u)\ r + I g z — p z I — uji max(g, p) 



(g - P ) 



2\2 



CO? 



3ujfq 2 p 2 



In 



ui + \q + p\ , (q 

ui + \q-p\ 



+ 



p 



2\2 



Qujfq 2 p 2 



uf uf + \q + p\ 2 - uJi\q + p\ 
+ \q - P\ 2 - ^l\q ~ P\ 



3gp ujf + \q-p\ 3 



I W- 
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dQ 



\q-p\ 2 +lo 2 + 2M 2 D » v y p ' q 

(f-p 2 ) 2 ln (g + p) 2 
2Ml+u 2 (q-p) 2 



2tt 7T 

= — + 



qp 2q 2 p 2 

{(q - p) 2 + 2M 2 D + oo 2 }{(2M 2 D + u 2 ) 2 + (g + p) 2 co 2 } (q + p) 2 + 2M 2 D + u? 



2M 2 D 



2M 2 (2M 2 +oj 2 ) 
ln 



(q-p) 2 + 2M 2 ) +^ 



{{q + p) 2 + uj 2 }{(q - p) 2 + uj 2 } (g + p) 2 + w 2 



(g — p) 2 + to 2 
0$tv(A±^A*Y) 



2n 7T 



D 

(q 2 -p 2 ) 2 ln (q+p) 2 



2M 2 D +lj 2 (q- p) 2 



qp 2q 2 p 2 

{(q + p) 2 + 2M 2 D + u; 2 }{(2M 2 + i^ 2 ) 2 + (q - p) 2 u 2 } ^ (q + p) 2 + 2M 2 + u? 



2Af£(2M£ + 



{q-p) 2 + 2M 2 D +oj 2 



{(g + P) 2 + ^ 2 }{(g ~ P) 2 + ln (g + P) 2 + ^ 2 
2M 2 n (q - p) 2 + to 2 



dQ 



\q — p} 2 + lo 2 



0( 3 Jtr(^ 7 %V)=vr 



\q-p\ 2 + L0 2 ■'■ ; " 



2 (g-p) 2 +o; 2 i (g + p) 2 +u; 2 

1 In 

qp 2q 2 p 2 (q — p) 2 + uj 2 

2 (q- pf + uj 2 (q+p) 2 +L0 2 

^ — O — O ^ 

gP 



2q 2 p 2 (q — p) 2 + a; 2 

(B-4) 



Here we list the products of the color matrices in the SDE. The completeness 
relation leads to 

8 1 1 

J~] (T ) a 'a{T ) Vb = -5 a >b5 a b> ~ T<W<W ■ (B-5) 
A=l 1 D 

By using the above formula the products of the color matrices in the SDEs (3-6) and 
(3-7) can be written 

5~ 



£ ( TA )aa'$ a b (T A )b'b = g^oft + S a3 5 b3 



A=l 



1-1 

^2(T )aa'5a'38b'3(T ) V b = ^ab + q $0,3^3 
A=l 1 6 



The product in the SDEs (3-8) and (3-9) can be written 



£(r A We^' 3 (T A M' = -^ 

A=l 6 



(B-6) 
(B-7) 
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To obtain the effective potential we need to take the trace and the determinant 
over the color, flavor, spinor and Nambu-Gorkov indices. The determinant in the 
first term of Eq. (3-1) is given by 



DetiiSpip)- 1 } = {F(p, B 1 ,A)} 8 {F(p, B%, 0)} 4 . 
The trace in the second term becomes 



(B- 



Tr{4 0) (p)- 1 ^(p)}=8 



F(p, Bi, A) 



2(p 2 + - p 2 )(pl - p 2 



-(pI-P 2 ){Bi(-p)} 2 - (p 2 - ~ p 2 ){Bi(p)} 2 
- ((Po) 2 - (p + ^ \A-\ 2 - ((p ) 2 - (p - p) 2 ) \A+\ 2 
(Po + A 1 ) 2 - P 2 



+ 



(Po - p) 2 - P 2 



(po+/U )2_p2_ {j B 3(p)} 2 (po _^2_p2_ {jB3( _ p)} 2 

(B-9) 



References 

1) G. E. Brown and M. Rho, hep-ph/0103102; Phys. Rep. 269(1996),333, hep-ph/9504250 

F. Wilczek, hep-ph/0003183 

S. P. Klevansky, hep-ph/9810399 
R. D. Pisarski, hep-ph/9503330 

T. Hatsuda and T. Kunihiro, Phys. Rep. 247(1994),221, hep-ph/9401310. 

2) B. C. Barrois, Nucl. Phys. B 129(1977),390. 

3) D. Bailin and A. Love, Phys. Rep. 107(1984), 325. 

4) M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B 422, (1998), 247, hep-ph/9711395 
R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Phys. Rev. Lett. 81(1998), 53,hep- 
ph/9711396. 

5) K. Rajagopal and F. Wilczek, hep-ph/0011333 
D. K. Hong, hep-ph/0101025 

M. Alford, hep-ph/0003185 

T. Schafer, Int. J. Mod. Phys. B 15(2001), 1474, nucl-th/9911017. 

6) D. K. Hong, V. A. Miransky, I. A. Shovkovy and L. C. Wijewardhana, Phys. Rev. D 61 
(2000),056001, [ibid. D 62, (2000),059903,Erratum], hep-ph/9906478. 

7) T. Schafer and F. Wilczek, Phys. Rev. D 60(1999), 114033, hep-ph/9906512. 

8) C. D. Roberts and S. M. Schmidt, nucl-th/9903075 

R. D. Pisarski and D. H. Rischke, Phys. Rev. D 61(2000),051501, nucl-th/9907041 

V. A. Miransky, I. A. Shovkovy and L. C. Wijewardhana, Phys. Lett. B 468(1999),270,hep- 

ph/9908212 

I. A. Shovkovy and L. C. Wijewardhana, Phys. Lett. B 470(1999), 189, hep-ph/9910225 
N. Evans, J. Hormuzdiar, S. D. Hsu and M. Schwetz, Nucl. Phys. B 581(2000), 391, hep- 
ph/9910313 

R. D. Pisarski and D. H. Rischke, Phys. Rev. D 61(2000),074017, nucl-th/9910056 
M. Matsuzaki, Phys. Rev. D 62(2000),017501, hep-ph/9910541 
K. Rajagopal and E. Shuster, Phys. Rev. D 62(2000),085007, hep-ph/0004074 
H. Abuki, T. Hatsuda and K. Itakura, hep-ph/0109013. 

9) J. C. Bloch, C. D. Roberts and S. M. Schmidt, Phys. Rev. C 60(1999),065208, nucl- 
th/9907086. 

10) J. Berges and K. Rajagopal, Nucl. Phys. B 538(1999), 215, hep-ph/9804233 
K. Langfeld and M. Rho, Nucl. Phys. A 660(1999),475, hep-ph/9811227 

G. W. Carter and D. Diakonov, Phys. Rev. D 60(1999),016004, hep-ph/9812445 



Phase Transition in Dense QCD with the Schwinger-Dyson Equation 



37 



A. Chodos, F. Cooper, W. Mao, H. Minakata and A. Singh, Phys. Rev. D 61(2000),045011, 
hep-ph/9909296 

H. Mishra and J. C. Parikh, Nucl. Phys. A 679(2001),597, hep-ph/0003019 

B. Vanderheyden and A. D. Jackson, Phys. Rev. D 62(2000),094010, hep-ph/0003150 
Y. Kim and M. Rho, nucl-th/0004054 

R. Nebauer and J. Aichelin, hep-ph/0101289. 

11) R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Ann. of Phys. 280(2000),35, hep- 
ph/9904353. 

12) D. T. Son, Phys. Rev. D 59(1999),094019, hep-ph/9812287. 

13) M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B 537(1999),443, hep-ph/9804403. 

14) M. Alford, J. Berges and K. Rajagopal, Nucl. Phys. B 558(1999),219, hep-ph/9903502 
T. Schafer and F. Wilczek, Phys. Rev. D 60(1999),074014, hep-ph/9903503. 

15) R. D. Pisarski and D. H. Rischke, Phys. Rev. D 60(1999),094013, nucl-th/9903023. 

16) N. Evans, S. D. Hsu and M. Schwetz, Nucl. Phys. B 551(1999), 275, hep-ph/9808444 Phys. 
Lett. B 449(1999), 281, hep-ph/9810514. 

17) T. Schafer and F. Wilczek, Phys. Lett. B 450(1999), 325, hep-ph/9810509. 

18) R. D. Pisarski and D. H. Rischke, Phys. Rev. Lett. 83(1999),37, nucl-th/9811104. 

19) V. A. Miransky, Dynamical symmetry breaking in quantum field theories (Singapore, Sin- 
gapore, World Scientific, (1993), p. 533). 

T. Kugo, "Basic concepts in dynamical symmetry breaking and bound state problems," 
KUNS-1086 Lecture delivered at 1991 Nagoya Spring School on Dynamical Symmetry 
Breaking, Nakatsugawa, Japan, Apr 23-27, 1991. 

20) T. Kugo and M. G. Mitchard, Phys. Lett. B 282(1992), 162. 

21) K. Higashijima, Phys. Rev. D 29(1984), 1228 
V. A. Miransky, Nuovo Cim. A 90(1985), 149. 

22) K. Aoki, M. Bando, T. Kugo, M. G. Mitchard and H. Nakatani, Prog. Theor. Phys. 
84(1990),683. 

23) J. M. Cornwall, R. Jackiw and E. Tomboulis, Phys. Rev. D 10(1974),2428 . 

24) Y. Taniguchi and Y. Yoshida, Phys. Rev. D 55(1997), 2283, hep-ph/9512375. 

25) M. Harada and A. Shibata, Phys. Rev. D 59(1999),014010, hep-ph/9807408. 

26) H. Pagels and S. Stokar, Phys. Rev. D 20(1979), 2947. 

27) J. Gasser and H. Leutwyler, Ann. of Phys. 158(1984), 142. 

28) J. Gasser and H. Leutwyler, Phys. Rep. 87(1982), 77. 

29) A. Barducci, R. Casalbuoni, S. De Curtis, R. Gatto and G. Pettini, Phys. Rev. D 
41(1990), 1610. 



